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(54) Procede et echographe pour la determination de la viscoslte et du gradient de pression dans 
un vaisseau sanguin 



(57) Proc6d6 de determination de la viscosity et du 
gradient de pression dans un vaisseau sanguin, com- 
prenant I'acquisition de n>2 valours de vitesses du sang 
en correspondance avec un meme nombre n de rayons 
du vaisseau sanguin, determines le long d'un diametre 
situe k une position axiale donn^e, formation d'un vec- 
teur de vitesses du sang au moyen de ces n valeurs de 
vitesses du sang, et evaluation de ladite viscosite et du 
gradient de pression ^ partir d'une transformation dudit 
vecteur de vitesses du sang avec tormation d'une rela- 



tion lindaire liant directement un vecteur des debits (y) 
au vecteur des derivSes de vitesses (h) factorise par la 
viscosite et au vecteur de gradient de pression (S), 
et evaluation simultanee des deux valeurs a determiner 
de viscoslte (p.) et de gradient de pression (5) k partir de 
cette equation directe. Les vitesses sont favorablement 
determinees par correlation de signaux ultrasonores 
foumis par un echographe. 

Application : Imagerie medicale, echographie me- 
dicals 
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Description 

[0001] L' invent ion concerne un proc6d6 de determination de la viscosity et du gradient de pression dans un valsseau 
sanguin. L'invention trouve son application dans le domaine de fabrication des appareils ^hographiques d ultrasons, 

5 appel^s ^hographes, pour I'aide au diagnoslic medical. 

[0002] Dans le domaine medical, lorsqu'un patient a un accident vasculaire cerebral ou cardiaque, il est n6cessalre 
de contrCler la viscosity du sang, car ce paramfetre joue un r6le important dans le diagnostic de I'accident vasculaire 
et dans la definition d'un traitement ult6rieur appropri6, par exemple un Iraitement k base de medicaments fluidifiants. 
Actuellement, pour connaTtre la viscosity du sang, il est n^cessaire de rdaliser des mesures sur des ^chantillons san- 

10 guins obtenus par des prdldvements sur le patient. II en rdsulte des inconvdniants importants. Le rdsultat des mesures 
n'est pas imm^diat, las prelevements eux-memes sont invasifs done potentiellement contaminants. Un autre inconve- 
nient est que les prei&vements doivent etre r^pet^s periodiquement pour controler I'etfet du traitement ce qui risque 
d'endommager les vaisseaux du patient. On ne connaU pas actuellement d'alternative k cette m^thode de determination 
de la viscoslte du sang. Par moyen invasif on entend un preievement du sang qui necessite par exemple Tintroductlon 

IS d'une aiguille dans le vaisseau sanguin pour emplir une seringue avec un echantillon de sang. 

[0003] II est ddj^ connu du brevet EP 0 430 374 (US 5 1 03 826. Bonnefous, Apr. 1 4, 1 992), un appareil echographique 
k ultrasons muni d*un systdme de traitement de signaux ultrasonores pour determiner des profile de la Vitesse du sang 
dans une artere. L'appareil 6chographique comprend un transducleur emettant et recevant des signaux ultrasonores 
d'une region d'artere d'un patient, et un etage transmetteur connects a la sonde. Le syst^me de traitement estime la 

20 Vitesse du sang dans fartere, en fonction du temps et de la profondeur d'exploration, par un proc6de de correlation 
de signaux successifs renvoyes par le transducteur. 

[0004] Un but de l'invention est de foumir un procede apte k etre mis en oeuvre, sans utiliser de moyens invasifs, 
pour permettre la determination de la viscosite du sang s'ecoulant dans une artere sous Taction de la pulsation car- 
diaque qui cree un gradient de pression, cette operation etant etfectue avec une precision egale ou superieure k celle 
qui etait atteinte par la methode connue incluant {'analyse de prelevements. 
[0005] Ce but est atteint par un procede selon la revendication 1 . 

[0006] Ce procede presente des avantages importants du fait qu'il est aussi precis que le procede par prelevements, 
qu'i! est plus rapide et plus simple k mettre en oeuvre. et qu'il peut etre repete un grand nombre de fois sans endom- 
mager les vaisseaux sanguins du patient, puisqu'il n'est pas invasif. Un autre avantage important reside dans le fait 
30 qu'aucun echantilbn de sang n'etant preieve ou traite, tout probieme de contamination du patient ou du personnel de 

laboratoire est evite. 

[0007] Un appareil echographique k ultrasons selon la revendication 10 est muni de moyens pour mettre en oeuvre 

ce procede. 

[0008] Le procede est decrit ci-apres en detail, en reference aux figures schematlques annexees dont : 

35 

la FIG. 1 qui represente sous fonme de blocs fonctionnels. un dtagramme de deroulement des etapes du procede 
de determination de la viscosite du sang dans une artere et du gradient de pression du k la pulsation cardiaque ; 
la FIG.2A qui represente un axe d'exploration echographique (Z) vis-^-vis de I'axe longitudinal (X) d'un vaisseau 
sanguin ; la FIG.2B qui illustre le decalage temporel entre deux tlrs echographiques et la FIG.2C qui illustre une 
40 fenetre de resolution; 

les FIGS.3A et 3B qui illustrent des premier et second signaux temporels et la FIG.3C qui illustre un signal de 
correlation ; 

la FIG. 4 qui illustre un signal de correlation en 1 bit; 

la FIG.5 qui illustre un profil de vrtesse le long d'un axe de tir echographique Z ; 
4S la FIG.6 qui represente sous forme de blocs fonctionnels un diagramme de deroulement du procede de transfor- 
mation d'un vecteur de vitesses pour foumir les valeurs de viscosite et gradient de pression dans le vaisseau 
sanguin. 

[0009] L'invention concerne un appareil d*echographie ultrasonore muni d'un sysieme de traitement des signaux 
50 ultrasonores pour la determination des parametres physblogiques du sang s'ecoulant dans une artere examinee in 
vivo dans une region tissulaire, et en particulier pour la determination de la viscosite du sang de maniere non invasive 
et en temps reel, et pour la determination du gradient de pression sanguine dans cette artere. 
[0010] En reference k la FIG.1 et aux FIGS.2A, B, C, cet appareil comprend au moins un transducteur ultrasonore 
100 associe k un etage d emission 200 et k un etage de traitement numerique 300. Le transducteur 100 est applique 
ss k la region tissulaire 11 pour I'examen de I'artere 10 ayant des parois 12, 13. Le faisceau emis par le transducteur 
ultrasonore est de preference large bande. 

[0011] L'etage d'emission 200 comprend un sequenceur compose d'un oscillateur et d'un divlseur de frequence qui 
commande, avec une periode de recurrence X un generateur dont les signaux eiectriques d'excitation sont envoyes 
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vers le transducteur 100 qui les convertit en signaux temporels ultrasonores, qui sont emis avec la pericKle de recur- 
rence T, selon un axe de tir Z. Dans ce systeme d' emission, chaque signal tempore! ultrasonore a un rang d'6mission 
dans le temps note k, k+1 , ...de periode T L'6tage d'6mission 200 reli6 au transducteur 100 assure ainsi la formation 
d'un faisceau 50 ultrasonore pour I'exploration de la region tissulaire 11 incluant Tart^re 10. Dans cet examen: on 

s suppose que I'^coulement sanguin se fait selon Taxe longitudinal X de I'art^re 10, at que Taxe de tir Z du faisceau 
ultrasonore ^ partir du transducteur 100 fait un angle 0 avec I'axe X. Les profondeurs d'exploratlon sont discretisees 
sur I'axe de tir Z, en des points not^s Rq, R^.... R|sj rep6r6s ^ partir d'un point 0. de coordonn^e spatiale sur Z not6e 
r^, situe ^ I'intersection des axes X et Z ^ I'lnterieur de I'artere, jusqu'a un point de coordonn^e spatiale sur Z not§e 
Tfg et situ6 sur une parol par exemple 12 de I'artere 10. 

10 [001 2] En reference ^ la Fl G. 1 , un 6tage sdparateur 1 01 est insdrd entre les dtages 200, 300 et le transducteur 1 00 
pour eviter I'aveuglement des etages 200 et 300 par les signaux d'6mission. L'etage de traltement 300 comprend en 
sortie de I'etage separateur 101 un amplificateur haute frequence incluant une compensatbn de gain en fonction de 
la profondeur d'exploratlon dans la region tissulaire examinee. L'etage de traitement 300 comprend aussi un eliminateur 
301 d'^chos fixes qui foumit, k partir d'un signal S^{X,r) provenant de l'etage separateur 101. un signal echographique 

IS tempore! 6^{rA) d^barrass^ des composantes fixes dues notamment aux reflexions sp^ulaires sur les parois 12, 13 
de I'artdre 10 examinee. 

[0013] Le signal 6^{r,\) issu de I'ellminateur d'echos fixes 301 est traite de preference selon une methode de corre- 
lation decrite dans le brevet EP 0 225 667 (US 4 803 390, Bonnetous). Cette methode est brievement rappelee ci-apr^s. 
[0014] Selon cette methode, le systeme de traitement 300 determine la Vitesse d'6coulement de structures du sang, 

20 par example des amas de globules rouges, appel^es cibles, notees TG (en anglais : Target), a I'instant t ou ces cibles 
coincident avec I'axe de tir Z et en des points discretises sur I'axe de tir Z appartanant ^ rinterleur de I'artere. II est 
suffisant d'estimer les vitesses entre le centre 0 et le point R|m sur la parol 12 car les vttesses entre le centre 0 et le 
point R'jg k I'intersection de I'axe Z et de I'autre parol 13 de Tartfere sont suppos6es sym6triques. Les profondeurs 
d'exploratlon le long de I'axe Z sont donnees par des coordonnees spatiales r Le temps ou un echo est renvoye par 

25 une cible situ^e k une profondeur r est note t. En particulier, le systeme de traitement 300 traite les signaux echogra- 
phlques temporels de rang k notes d|^ (r,t) formes des Intensltes des echos renvoyes par les cibles en chaque point 
discretise de I'axe Z de profondeur ro k r^ dans I'artere 10, et foumit la vrtesse d'ecoulement du sang en ces points. 
Les FIGs. 3A et 38 representent de tels signaux temporels echographiques d|((r,t) et d|(^.-,(r,t) avec en ordonnee I'in- 
tenslte du signal echographique en unite arbitraire et en abscisse des coordonnees temporelles t en microsecondes 

30 qui correspondent aux coordonnees spatiales r le long de I'axe de tir Z. 

[0015] Les FIGS.2A et 2B iliustrent que rappareil echographique n'est pas sensible aux vitesses autres que celles 
qui se rapportent k I'axe de tir Z. Done, on ne considere ni les distances qui sont en realite parcourues par les cibles 
TG selon I'axe X, ni les vitesses de ces cibles selon cet axe X, mais seulement leurs projections sur I'axe Z. La vitesse 
d'une cible est done mesuree par sa composante V2(r,t] le long de I'axe Z en fonction du temps t et de la profondeur 

35 d'exploratlon r. Sur la FIG.2C, on a note ri^ et ri^^^ les profondeurs d'une cible TG renvoyant des echos dans les signaux 
de rang k et k+1 respectlvement. La composante selon I'axe longitudinal X se deduit de V^Crt) par la relatk}n : 



Vx(r,t) = V2(r,t)/cos0 (1) 

40 

[0016] Le systeme de traitement 300 met en oeuvre une technique de correlation de deux signaux temporels echo- 
graphiques successifs notes de rang k et k+1 . La technique de correlation demande, en reference k la FIG.2C: qu'un 
fort pourcentage de cibles present dans le faisceau ultrasonore lors d'un premier tir de rang k reste present dans le 
faisceau echographique tors d'un tir ulterieur de rang k+1 , de manlere k ce que ia correlation entre les deux signaux 
temporels successifs soit sufTisamment forte. C'est pourquoi, en reference k la FIG.2B, la correlation est effectuee 
dans une certaine fenetre de resolution notee RG de dimension w (en anglais: Range Gate) choisie pour que cette 
demande soit facilement satisfaite. Par exemple, si le sang a une vitesse de flux de I'ordre de 1 m/s, si la periode T 
entre deux excitations est egale k 100 ^is ce qui correspond k une frequence d'excltatlon de I'ordre de 10 kHz, le 
deplacement des cibles entre deux tirs successifs est de I'ordre de 0,1 mm. Ce deplacement est substantial iement 
petit devant les dimensions d'une cellule de resolution RG avantageusement choisies de I'ordre de 1 mm sebn I'axe 
X et 0,3 mm selon I'axe Z. Dans ce cas, environ 20 % de cibles quittent la cellule de r6solutbn pendant la periode T 
environ 20 % y entrent et un fort pourcentage de cibles reste present. 

[0017] L'idee basique du traitement mis en oeuvre par le systeme 300, pour la determination de la vitesse en un 
point discretise de I'axe Z, est qu'un signal tempore! ultrasonore cl|<+^(r,t) retrodiffuse par une cible en mouvement, 
sHuee k I'lnterieur de la fenetre de resolution RG, qui a ete emis par un tir echographique selon Z, de rang k+1, au 
temps t. est la repllque d'un signal tempore! anterleur emis par un tir echographique de rang k selon le mdme axe Z 
au temps t -T, retrodiffuse par la meme cible au temps t - T x, ou t est le temps supplementaire appeie decalage 
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tempore! (en anglais : Time Shift) mis par le signal tempore! ultdrleur pour parcourir !e trajet sonde-cible-sonde. du fait 
du d^placement de la cible durant I'intervalle T. 

[0018] Ainsi, au premier ordre, deux signaux temporels ^chograpliiques successifs retrodlffus^s par la meme cible 
sont \\6s par la relation : 

dk.i(t) = d,(t-T-T) (2) 

[0019] D'ou il r^ulte que !e d^calage tempore! x est donn6 par : 

X = 2V2(r.t)/c (3) 

ou c est la Vitesse du son dans le vide et par extension dans les tissus. 1^ lormule (3) permet de determiner ies 
composantes V2(r,t) des vrtesses des cibles en cliaque point discretis6 d'ordonn^es spatiales r en faisant une approxi- 
mation basee sur le fait que la vitesse axiale V^Cnt) est tres petite devant la vitesse du son c, et que le d^placement 
d'une cible est tr5s petit devant !a distance du transducteur k la cible durant une p6riode T 

[0020] Par consequent, !a vitesse du flux peut etre mesuree par un appareil ecliographique ayant un sysl^me de 
traitement 300 munis de moyens pour mesurer le decaiage tempore! x induit par les cibles entre deux tirs dchographi- 
ques. 

[0021] Au premier ordre. la fonction d'intercorr^lation entre les deux signaux temporels successils est de la forme : 

Qck.i(tb,u) = £^^M0 d,,,(t'^u)dt (4) 

ou le temps tQ est Ii6 ^ la profondeur d'exploration r par to= 2t/C: Wq est la largeur de la fenetre de resolution sur Taxe 
du temps, et u est un parametre de decaiage tempore! de la fonction d'intercorr^latlon (en anglais : cross-correlation). 
[0022] En introduisant I'equation (2) dans !'6quation (4). il vient : 



Ck,k.i(to.u) = C,k(to.u-x) (5) 

[0023] La fonction Ci^j^ (t^.u - x) est une fonction d'autocorrSlation. et, ce fait est maximale pour u = x . Dans une 
hypothese d'une situation ideale de vitesse axiale unifonme, le maximum de la fonction de correlation correspond au 
decaiage de temps x et par cons^uent permet une determination de la vitesse du sang pour le tir ^chographique au 
rang k selon la relation : 

40 Vz(r,t) = (c/2)(x^) (6) 

En reference k la FIG.1, le circuit 302 d'intercorrelatlon du syst6me de traitement 300 determine les vltesses V2(r,t) 
en cherchant pour quel parametre u la fonction (^j^^^ (Xq,u) est maximale. A cet effet, la fonction d'intercorrelation est 
6chantillonnee k un pas d'echantillonnage At, de fagon k obtenir 2m+1 valeurs de la fonction de correlation, m etant 
un scalaire. La valeur maximale de ces 2m+1 valeurs correspondant ^ u = Uq permet de mesurer x en utillsant I'dgalitd : 
T = Uo (7). 

[0024] Selon cette methode, il n'y a pas d'ambiguTte k la mesure des vitesses parce qu'il n'existe qu'un seul maximum 
k cette fonction de correlation. 

[0025] La methode de decaiage temporel derive du deplacement reel des cibles entre deux tirs successifs dans le 
temps selon le meme axe Z. C'est pourquoi la fenetre de resolution doit etre habilement choisie. La FIG.3A illustre 
une premiere fenetre de resolution de dimension Wq = 2 )iS appliquee k un signal temporel d|((r,t) et ia FIG.4B illustre 
une seconde fenetre de resolution de dimension = 4 appliquee au signal ulterieur d|(^i(r,t). La FiG.4C Illustre la 
correlation C des signaux temporels d|^(r,1) et d|(^^(r,t) qui conduit k la determination d'un pic Ci^^x^e correlation pour 
x de I'ordre de 0,3 |j,s ce qui permet de calculer la vitesse V^trt) par la relation (6). 

[0026] Afin de s'affranchir des en'eurs inherentes k rechantillonnage dans la determination du maximum de la fonc- 
tion de correlatbn, le systeme de traitement 300 comprend favorablement un circuit 303 de multiplexage-interpolation 
fournissant k partir des valeurs des fonctions de correlation une estimee plus precise de valeur du pic de correlation 
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et done de ia Vitesse correspondante. Dans ce type de traitement, la correlation entre signaux est une correlation dite 
«1 bit » en ce sens que les signaux yk+i(i',t) et yk(r,t) utilises pr6c6demment sont reduits ^ leur signe. On sail que, dans 
ce cas, le pic C|^;^ de la fonction de correlation est de forme triangulaire tsoc^le. Comme illustr^ par la FIQ.4, la 
connaissance de cette forme permet de n'eff ectuer le calcul de correlation qu'en des points discrets des signaux tem- 
s porels, et, une fois determines un certain nombre de points, permet k partir du point le plus haut et de ses deux volsins 
de reconstrtuer compl^tement, par interpolation lineaire. le pic de correlation Cf^^^ et done de determiner avec precision 
remplacement de Uq = x, d'ou Ton peut d6duire VzCr.t). 

[0027] En reference h la FIG. 1 , les valeurs trouv^es pour la vitesse V(t,z) sont stockees dans une memoire 304 de 
stockage en vue de traitements ulterieurs. Comme illustre par la FIG.5, en reportant les valeurs de vitesse V2(r.t) en 
10 cm/s en ordonnde et les profondeurs r en mm en abscisse, on obtient un profil des vitesses sur I'axe Z dans I'art^re 
au temps t. De meme que I'axe Z est discretise en N points d'une part de O et en N points tfautre part de O, le profil 
des vitesses est discretise avec le meme nombre de points. 

[0028] En reference k la FIG.1, k I'aide de moyens de visualisation connus, il est possible d'afficher sur un ecran 
305 une image en couieur de reeoulement du sang qui est derivee de revolution au cours du temps du profil de vitesse 

IS dans le vaisseau. Le codage couieur de la vitesse dans cette image est notee CVI (de I'anglais ; Color Velocity Imaging). 
[0029] Un avantage important de la methode echographique de correlation en 1 bit est que les profils des vitesses 
sont obtenus en temps r6el. Cependant toute autre methode apte a fournir des donnees de vitesse du sang dans un 
vaisseau avec une precision suffisante peut §tre utilisee pour mener d bien revaluation de la viscosity et du gradient 
de pression comme decrit ci-apres. 

20 [0030] Les profils de vitesse du sang, k la m§me localisation Z, varient dans le temps. Les caracteristiques d'un 
fluide telles que la viscosity, la compressibilite, Tisotropie ; la nature du flux telle que laminaire ou turbulente ; les 
proprietes du tube qui contient le fluide, telles que la rigidite, reiasticrtS, ou la viscoelasticite ; et la forme du gradient 
de pression telle que large-bande ou bande etroite, sont des param6tres qui dictent la forme des profils de vitesse. 
[0031] On suppose ci-apres que le sang dans une artere est assimilable a un fluide Newtonien ayant un flux laminaire 

25 s'ecoulant dans un tube rigide. Comme illustre par la FIG.5, sous un gradient de pression note dP/ dx, la vitesse V(r, 
t) suit une bi de Navier-Stokes, exprlmee par la tormule connue (8) : 

dP ,dV. ,(fv. 1 dV 
ou p est la densite du fluide et ou ^ est la viseosite 

[0032] La loi de Navler-Stokes, qui lie la viseosite et le gradient de pression k la vitesse, ne permet cependant pas 
de calculer directement ia viseosite du fluide et le gradient de pression en fonction de la vitesse du fluide. Le calcul de 

SB la viseosite du fluide par inversion de I'dquation de Navier-Stokes conduit k des expressions analytiques complexes. 
[0033] Selon ('invention, ces difficultes sont evitees, et la viseosite et le gradient de pression sont determines en 
temps reel par un procede de traitement de donnees incluant une inverskxi de la relation (8) qui consiste k estimer la 
viseosite et le gradient de pression a partir des donnees fournies par des profits de vitesse prealablement determines 
de preference comme decrit pr6c6demment 

40 [0034] Sous I'hypothese d'un flux laminaire, requation de Navier-Stokes gouveme les variations de la vitesse V(r,t) 
pour r et t donnes. Pour les points discretises Rj de Rq ^ Rn> coordonnees spatiales rj de rg k r^ sur I'axe Z, 
I'integration de {'equation (8) de Navier-Stokes entre r = r^ et r = r^ foumit {'expression (9) sulvante, ia derivee radiale 
de Vx etant nulle : 

ncj dt " dx dr 

ou rj est la profondeur d 'exploration 
so [0035] Chaque membre de requation (10) est une fonction du temps t et du rayon r^. Pour chaque echantillon de 
temps, la relation peut dtre ecrite sous une forme vectorielle, chaque vecteur etant forme par les valeurs des differentes 
quantites k differentes profondeurs d'exploratlon rj. D'ou : 
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(10) 



IS 



Liquation du probl^me direct peut done etre mise sous la torme : 

y = jih + 6 + b 



(11) 



20 



oCi b est un bruit blanc gaussien qui modeiise les erreurs de mesure. En accord avec I'equatlon (11), y est un vecteur 
gaussien dont la moyenne est egale d : 



)ih + 6 



(12) 



[0036] La phase essentielle du proc^dd constste ^ Inversercette Equation qui mod6lise le probldme direct (11). La 
solution du probleme inverse est obtenue de preference au moyen d'un esttnnateur de maximunn de vraisemblance 
2S connu (Maximum Likelihood Estimator) utilise pour estimer k la fois la viscosity [i et le gradient de pression 5. 

[0037] On cherche les valeurs p. et 5 qui maximisent la vraisemblance de y connalssant ^ et 5 donn^es par une 
formule (1 3) : 



30 



(Ji,S) - argmax^^^[p(y|/^,(y)] (13) 



35 



ou p(ylp„ 5) est la probabilite de y connaissant |i, 6, et oij |i et 5 sent des estimateurs de |i et 5 appelees estimateurs 
du maximum de vraisemblance. Ceci est equivalent a rechercher le minimum d'une fonclion de cout notes W(jj.,6): 

W(ji,6) = (y - - 5)"^ Kl (y - ^Ai - 6) (14) 



ou K'1„ est une matrice diagonale dont les coefficients dependent de la confiance que Ton accorde aux valeurs de y, 
40 ce qui conduit aux expressions sulvantes de m et 5 : 



45 



A 



r '^ r -1 
y ^nh- — — 

1 Vnl 

h k„ n — j 

1 %1 



(16a) 



so 



^= ~^ 

/I Kn1 



(17a) 



55 [0038] En premiere approximation, on suppose que le bruit est distribu^ de fagon identique, de mani^re k ce que la 
matrice K„ soil proportionnelle k la matrice identity. Ceci produit une premiere approximation de la solution : 
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40 



45 



SO 



y h- — f— 
0(0) = 1-V 



6(0)= ^V-A^^^ (17b) 



70 /? 1 

[0039] Selon I'invention, on propose un proced6 de determination de la viscosity \i fonde sur cette methode de 
resolution du problems inverse, a partir du problems direct constitud par la relation (11) et des profils de vitesse de- 
termines de preference par echographte. ou par toute autre methode donnant des resultats equivalents, en utilisant 

IS cet estimateur de maximum de vraisemblance selon des etapes decrites cl-apres. 

[0040] En reference aux FIG. 1 et 6, un processeur 400 est associe au systems de traitement 300 de Tappareil 
echographlque ou d'un systdme fournissant les vitesses. En 6chographie, le systeme de traitement 300 est «materiel 
>> au sens HARDWARE utilise par I'homme du metier, et le processeur 400 est programme au sens SOFTWARE 
egalement connu. Le processeur 400 peut faire partie de rechographe ou bien d'une station de travail coupiee k 

20 rechographe. 1.8 programmatlon du processeur 400 definit les differentes etapes necessaires pour determiner la vis- 
costte \i et le gradient de pression 5 sebn le precede de {'Invention illustre par les blocs fonctionnels de la FIG.6. Ce 
procede comprend des etapes de : 

acquisition 410 des profils de vitesses Vx(r,t) k partir des profils de vitesse V2(r,t) foumies par le systeme de 
2S traitement 300, en utilisant la fomnule (1 ) ; 

echantillonnage 420 d'un profil de vitesse discretise donn6, en utilisant un nombre n de points Ri Si Rn parmi les 
N points de discretisation a R|^, tel que 2 < n < N ; les n points de segmentation definissent des rayons r^ a r„ 
du tube assimile au vaisseau sanguin, notes de manlere generate rj; 

determination 430 d'un vecteur de vitesses constitue des vitesses aux points de segmentations R^ k R^, ces 
30 vitesses etant notees \^x(h) ^ Vx(rp) et d'une maniere generale Vx(rj). 

evaluation 440 des debits du fluide dans le tube, c'estnd-dire du sang dans le vaisseau, dans les sections de tube, 
definies par les rayons rj de r^ a r,,, ce qui fournit respecllvement les debits notes Q(rj) de Q(r^) a Q(rn) ; 
calcul 450 du vecteur y selon la fomriulation (10) en etfectuant des calculs des quantites 



nfj dt 



(20) 



pour Tj de ri ^ r^ ; 



Ces calculs (20a) comprennent des sous-etapes de : 
determination du debit Q(rj) calcuie selon la formulation sulvante : 



Q(rj)-y2;rriVV(r,) (20a), 



determination des derivees temporeiles des debits Q(rj) selon : 



dO(r.) 

pour r^ de r^ k r„ (20b), 



ss et multiplication des derivees temporeiles des debits par les quantites respectives avec q de ^ : 
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pMr^ (20c) 

ou p est la masse volumique du fluide suppos^e environ 6gale k 1 ; 

calcul 460 du vecteur h selon la fomiulation (10) en effectuant des calculs des quantrt^s : 



! (30) 



10 'Z 

qui sont les d6riv6es axiales des vitesses pour chaque point R, de Ri ^ d6finissant les rayons fj ; 

calcul 470 d'un produit matriciel : K'^ M. ou la matrice M est donn6e par la formule : M= / - p.rt - 5, et ou M"^ est 

la matrice transpose de M. 

calcul 480 de la viscosit6 \i selon les formulations (16) 
calcul 490 du gradient de pression 5 = -^elon (17) 

[0041] Pour etfectuer les stapes 410 ^ 490, le processeur est programme par les fonmules de d6riv6es temporelles 
de d6bit, de d^riv^es axiales de Vitesse, et de calcul matriciel. Le syst^me de traitement 300 fournlt automatiquement 
en temps r6el les profils de vitesse et le processeur 400 fournit automatiquement subsdquemment la viscoslte et le 
gradient de pression dans un delai total extremement court c'est-^-dire en temps r6el. 

[0042] Grace au precede selon Tinvention. pour resoudre I'equation (8) qui a deux inconnues \i et 5, on procede 
done par une inversion du probl6me direct qui utilise n points discretises d'un profil de vitesse pour resoudre en realite 
n Equations k deux inconnues, avec n > 2. Dans le cas de la determination de la viscosite du sang, il est interessant 
de selectionner un nombre n > 2 de points de discretisation pour disposer de plus de 2 equations independantes a 
deux inconnues pour la determination des 2 inconnues. En effet, si Ton se reduit a n = 2, on peut se trouver occasion- 
nellement devant une indetermination qui ne permet pas de determiner les deux inconnues. Dans ces conditions les 
moyens de calculs sont bloqu6s. Alors que si Ton choisit favorablement un nombre n>2 de points de segmentation, 11 
se trouve toujours un nombre sufRsant d'equations pour determiner les deux inconnues p. et 5. 

30 



Revendlcatlons 

1. Proc6d6 de determination de la viscosite et du gradient de pression dans un vaisseau sanguin. comprenant des 
etapes de : 

acquisition d'une plurality n>2 de valeurs de vitesses du sang en correspondance avec un m§me nombre n 
de rayons du vaisseau sanguin, determines le long d'un diametre situe k une position axlale donnee, et k 
differents instants, 

formation d'un vecteur de vitesses du sang au moyen de ladite pluralite de valeurs de vitesses du sang, 
evaluation de ladite viscosite et du gradient de pression k partir d'une transfomnatlon dudit vecteur de vitesses 
du sang. 

2. Precede selon la revendication 1 , comprenant, pour effectuer la transfomnatlon dudit vecteur de vitesses du sang 
etapes de : 

fomnation d'un vecteur des debits (y) constitu6 d'une mSme pluralite n>2 de valeurs de debits du sang k travers 
des sections du vaisseau sanguin en correspondance avec les n valeurs de rayons, 
formation d'un vecteur (h) des derivees des vitesses constitue d'un quotient de la meme pluralite n>2 de 
derlvees radiales des vitesses par les rayons correspondants, ce vecteur etant factorise par la viscosite (p.) k 
determiner, 

formation d'un vecteur de gradient de pression (5) constitue du produit du gradient de pression k determiner 
par une matrice identite, 

fomnation d'une relation lineaire liant directement le vecteur des debits (y) au vecteur des derivees de vitesses 
(h) factorise par la viscosite et au vecteur de gradient de pression (5), 

evaluation simultanee des deux valeurs k determiner de viscosite (^) el de gradient de pression (5) k partir 
de cette equation directe. 



4S 



so 



ss 
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3. Proc6d6 selon la revendication 2. comprenant revaluation simuttanee des deux valeurs k determiner de viscosity 
ip) et de gradient de pression (5) par une operation d'inversion de la relation lin^aire directe liant le vecteur des 
d^brts (y) au vecteur des d^riv^es de Vitesse (h) factorisS par la viscosity (^) et au vecteur gradient de pression (5). 

s 4. Procddd sebn la revendication 3, comprenant, pour effectuer I'inversion de la relation lin^aire directe liant le vecteur 

des debits (y) au vecteur des derivdes de vitesse (h) factorise par la viscosite (p.) et au vecteur gradient de pression 
($), une etape d'estimation du maximum de vraisemblance du vecteur des debits (y) connaissant la viscosite {p) 
et le gradient de pression (5). 

10 5. Procdde selon la revendication 4, comprenant. pour estimer ledit maximum de vraisemblance du vecteur des 
debits (y) connaissant la viscosite (^) et le gradient de pression (5) des etapes de : 

evaluation k partir de ladite relatbn lineaire directe d'une premiere matrice (M) telle que: 

IS 

M = y - p,h - 8 

et d'une seconde matrice transpose de la premiere matrice (M^), 
evaluation d'une troisieme matrice (K) qui est une matrice diagonale dont les coefficients dependent de la 
20 distribution du bruit, 

et evaluation dudit maximum de vraisemblance du vecteur des debits (y) connaissant la viscosite (p.) et le 
gradient de pressbn (5) en effectuant une recherche du minimum d'une fonction de coQt (W) qui est le produit 
de la premiere, troisieme et seconde matrices tel que : 

-IT 



6. Procede selon la revendication 5, ou la troisieme matrice diagonale (K) est une matrice identite diagonale, le bruit 
etant considere uniformement reparti. 

30 

7. Procede selon la revendication 2, ou le nombre des valeurs de vitesses est deux (n=2) et ou revaluation des deux 
valeurs k determiner de viscosite (ji) et gradient de pression (5) comprend des etapes de : 

lormation d'un couple de relations lineaires directes liant le vecteur des debits (y) au vecteur des derlv6es 
de Vitesse (h) factorise par la viscosite (^) et au gradient de pression (5). de la forme 

35 

y = ph + 5 



resolution du couple de relations lineaires directes, au moyen des deux valeurs de debits et de deux valeurs de 
40 derivees des vitesses, pour foumir la viscosite (p) et le gradient de pression (5). 

S. Procede selon Tune des revendications 1 k 7 comprenant pour I'acquisition de valeurs de vitesses des etapes de : 

acquisition de signaux temporelsfournis par un transducteur ultrasonore explorant un vaisseau sanguin selon 
^ un axe de tir echographique (Z), 

traitement des signaux temporels dans le domaine temporel pour fournir des signaux conreies temporellement 
et un profit de vitesses du sang discretise le long de la direction de Faxe de tir, 

estimation d'un vecteur de vitesses du sang a partir d'un nombre n > 2 de valeurs de vitesses seiectionnees 
sur le prof il des vitesses a dlfferentes profondeurs d'exploration (r) sur I'axe de tir (Z) dans le vaisseau sanguin. 

so 

9. Procede selon la revendication 8, ou retape de traitement des signaux echographiques temporels est une corre- 
lation temporelle en 1 bit. 

10. Appareil echographique pour la mesure et I'affichage de parametres physiologiques du sang s'ecouiant dans un 
ss vaisseau sanguin, incluant des moyens de mesure de la viscosite et du gradient de pression du sang, ces moyens 

comprenant : 

des moyens (200) pour I'acquisition de signaux temporels fournis par un transducteur (100) ultrasonore ex* 
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plorant un vaisseau sanguin selon un axe de tir ^chographtque (Z), 

des moyens (300) de traitement des signaux temporels dans le domaine temporel pour toumir des signaux 
corr^l^s temporellement et un profil de vitesses discr6tis6 le long de la direction de I'axe de tlr (Z). 
des nrx)yens (400) pour : 

restimation d'un vecteur de Vitesse du sang h partir d'un nombre n > 2 de valeurs de vitesses s^lectionn^es 
surie profil des vitesses ^d iff ^rentes prof ondeursd'exploration (r) surPaxedetir (Z) dans le vaisseau sanguin, 
et I'dvaluatton de la vlscosit6 du sang i\x) et du gradient de pression (5), par une transformation dudtt vecteur 
des vitesses en utillsant un proc^dS selon Pune des revendicatlons 2^9. 
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